Skeleton expansions for directed polymers in disordered media 
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Partial summations of perturbation expansions of the directed polymer in disordered media 
(DPRM) enables one to represent the latter as skeleton expansions in powers of the effective cou- 
pling constant A(£), which corresponds to the binding state between two replicas in the replica 
field theory of DPRM, and is equivalent to the binding state of a quantum particle in an external 
5-potential. The strong coupling phase is characterized by the exponential dependence of A(t) on t, 
A(t) ~ exp(p c t) with p c being the binding energy of the particle. For dimensions d > 2 the strong 
coupling phase exists for Aq > A c (d). We compute explicitly the mean-square displacement and 
the 2nd cumulant of the free energy to the lowest order in powers of effective coupling in d = 1. 
We argue that the elimination of the terms exp(p c t) in skeleton expansions demands an additional 
partial summation of skeleton series. 
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The behavior of a directed polymer in a disordered 
medium is a fundamental problem in statistical physics, 
which is related to different problems such as non equi- 
librium dynamics of surfaces, growth processes such as 
Eden model, growth by ballistic deposition, randomly 
stirred fluids (Burgers equation), dissipative transport in 
the driven-diffusion equation, and the behavior of flux 
lines in superconductors [l|-[3]- Despite enormous inter- 
est over the recent years (for reviews see 0-0) there 
is no basic understanding on the analytical description 
of the strong coupling behavior of the directed polymer. 
Moreover there is even a controversy on the existence 
of the upper critical dimension d c = 4 [3] -[Hj or infin- 
ity [li|-[2^. Very recently a connection between DPRM 
and random matrices was established in [§| . It was also 
found that DPRM possesses a multicritical behavior at 
the transition in high dimensions ll|-[icj. 

The distribution function of the free polymer end 
W(x\, t) the equation 



(1) 



where the random potential /i(x, t) is Gaussian dis- 
tributed and possesses the moments 

Mx,t)] =0, [/x(x,%(x',t')] = A 6 d (x-x.')6(t-t'). 

The Cole-Hopf transformation W(x, t) = 
cxp (^h(x.,t)j maps (UJ to the Kardar-Parisi-Zhang 
equation 0] 

^ = J v 2 h + ~(Vh) 2 + r ] (x,t), (2) 

where /i(x, t) is a single-valued function, which describes 
the height profile above a basal d-dimensional substrate 
x in the co-moving coordinate system, A is responsi- 
ble for the lateral growth, 7 is the surface tension, 
and the noise ?y(x, £) — (2^ / X)n(x,t) with the correla- 
tor (ri(x,t)rj(x' ,t')) = 2DS d (x - x')5(t - t') and D = 



2Ao7 2 /A 2 . The introduction of the new variable v 
results in the Burgers equation 0] 

^ = 7V 2 w fe + \viV k Vi + / fc (x, t) 



Vh 



(3) 



with the stirring force f (x, t) = Vr/(x, t). 

The behavior of the directed polymer is determined by 
two competing factors: i) Energy win from low valued 
sites of the random potential, which requires a transver- 
sal wandering of the polymer and ii) the lost of clastic 
energy due to polymer stretching. The critical exponents 
C and u>, which describe the behavior of the transversal 
displacement of the free polymer end and the fluctuation 
of the free energy 



[<x 2 {t) >] ~t 2C , 8F{t)~t u 



(4) 



lu = 1/3 11 . The 
1, which is due to 



are exactly known for d = 1: £ = 2/3. 
exponents fulfil the relation 2£ = u> + 
the Galilean invariance. 

The motivation for the present approach is the follow- 
ing. The renormalization group (RG) method permits 
to reorganize the bare perturbation series as expansions 
in powers of the effective coupling constant. In the case 
of DPRM it was shown that the one-loop RG is exact 
[13]- but the effective coupling constant increases 
under renormalization, and corresponds to the so-called 
run-away situation. The effective coupling constant pos- 
sesses a pole at a finite length, which restricts the renor- 
malization until this length. In this paper instead of us- 
ing RG we perform partial summations of perturbation 
expansions of DPRM and arrive at the skeleton expan- 
sions in powers of the effective coupling constant A(t), 
which corresponds to the exact four-vertex of the replica 
field theory of DPRM The effective coupling constant 
A(t) is well-defined for all t, and increases exponentially 
with t, A(t) ~ exp(p c i), with p c being the binding en- 
ergy in the 5-potential. The strong coupling phase man- 
ifests itself in the exponential increase of A(t) with t. 
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For dimensions d > 2 the strong coupling phase exists 
for A > A c (d). We compute explicitly the mean-square 
displacement and the 2nd cumulant of the free energy to 
the lowest order in effective coupling in d = 1. A funda- 
mental problem of elimination of the exponential terms 
exp(p c t) from skeleton expansions remains, however, un- 
solved. Relying on the problem of the perturbational 
study of the quantum particle in an external (5-potential 
we argue that the compensation of the exponential terms 
demands an additional partial summation of skeleton se- 
ries, which physically corresponds to the redefinition of 
the ground state. 

The free energy and its fluctuations can be obtained 
directly from the replica partition function [Z n ] , which is 
the partition function of n quantum particles with Hamil- 
tonian 

a=l a,b=l 

The perturbation series of [Z n ] in powers of Ao can be 
represented by graphs examples of which are shown in 
Fig. [TJ The transformation 

[Z n ] = ^^+^11 



The computation of the first two cumulants for polymer 
with fixed one end to the lowest order yields 

r A 1 .l-d/2 _ 

2 (2717)^/2 l-d/2 

A direct computation shows that for a polymer with both 
free ends the disorder is irrelevant. 

Let us now consider the connected part of the four ver- 
tex [Z 2 ] which is shown in Fig. [2] The graphs describe 
the interaction between two polymers (quantum particles 
for imaginary times). As in the two particle problem the 
center of mass and relative coordinates can be separated 
by using the appropriate transformation of coordinates. 
The part of [Z 2 ~\ c depending on relative coordinates coin- 
cides with the perturbation series of the Greens function 
of a polymer in an external (5-potential. This series deter- 

[z 2 ]=^H + T^ + Xll +... 

FIG. 2: The perturbation series of the four vertex in powers 
of A . 



FIG. 1: Examples of graphs associated with the perturbation 
series of the replica partition function. 



[Z n ] = [e nlnZ ] 



-nF/k B T 



dFp(F)e 



-nF/k B T 



shows that the inverse Laplace transform of the replica 
partition function [Z n ] yields the distribution function of 
the free energy p(F) [26[. The cumulants Ck of [Z n ], 
which are defined by the equation 



[Z n ] 



I c 

ex p(E 

Vs=l / 



read 



C\ = [hxZ], C 2 = [ln 2 Z] - [InZ] 2 , 

C 3 = [In 3 Z] - 3 [In 2 Z] [In Z\ + 2 [In Zf . 



Thus, the free energy and its fluctuations (in units of 
ksT) are given by cumulants of the replica partition func- 
tion. 

The transversal displacement of the free polymer end 
with one end fixed at the origin is computed to the lowest 
order in A as 



[< x 2 (t) >] = dr/t(l + 



1 



Ao 



4 - d (4ir-y) d / 2 



t l - d ' 2 + •••)• 



mines the renormalization of the bare coupling constant 
Ao. The straightforward summation of the perturbation 
series in d < 2 gives 



A(p) 



Ao 



1 - Ao(47r 7 )- d / 2 r(l - d/2)p-( 1 - d / 2 ) 



(5) 



where T(l — d/2)p ^ d / 2 ^ is the Laplace transform of 
the probability G o (0, t; 0) ~ t~ d/2 . The pole of © deter- 
mines the binding energy p c of the quantum particle in 
(5-potential. The inverse Laplace transform of ([5]) is well- 
defined for all times t. For large t in d = 1 we obtain 
A(t) as 



A(i) = ^l e * A o/47 + /l r 3/2 + ... 

27 \ TT 



(6) 



Eq. {5}, as it is, applies for d < 2 only. To extend the 
applicability of Eq. ([5] to dimensions d > 2 we replace 
the return probability t~ d l 2 by i~ d / 2 exp (— A 2 /i), where 
A is a microscopic cutoff. The effective coupling constant 
is given in this case by the expression 

A(p) = 



1 - 2A (47r 7 )- d / 2 



K, 



d/2-1 



(2X^p) 



(7) 

where K n (x) is the modified Bessel function of the sec- 
ond kind. For d < 2 Eq. ([7|) coincides with Eq. (0 in 
the limit A — + 0. The necessity of the cutoff in d = 2 
follows from the connection of DPRM to quantum parti- 
cle in a weak potential well, where in d — 2 the binding 
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energy depends both on the depth and the width of the 
potential well, and not as in d < 2 only on the prod- 
uct of depth and width. The analysis of the eigenvalue 
equation for d > 2 (the denominator of Eq. ([7])) yields 
that a binding state exists in dimensions, if Ao exceeds 
a threshold value Ag(d). Consequently, in accordance 
with Eq. ([6]) the effective coupling constant increases ex- 
ponentially with time for all dimensions, if Ao exceeds 
the threshold value Ag(d). The latter suggests that the 
strong coupling phase exists for d > 2, if the condition 
A > Ag(d) applies. 

We now will perform the partial summation of pertur- 
bation series, which take into account insertions associ- 
ated with the [Z 2 ] -vertex shown in Fig. [2j and results in 
skeleton graphs, where the bare interaction is replaced by 
the effective coupling constant A (£2 — ti)- The possibil- 
ity for such a resummation is based on the fact that the 
numerical prefactors of the graphs under resummation 
are the same. Figs. [30 give examples of skeleton graphs 
contributing to the mean-square displacement and the 
cumulants C 2 . 



[<X 2 (t)>] = 



Carrying out the integrations over space variables results 



I 




FIG. 3: Examples of skeleton graphs contributing to the 
transversal displacement. 



c 2 = 



FIG. 4: Examples of skeleton graphs contributing to the 2nd 
cumulant of the free energy. 

The expressions for the mean-square displacement and 
the cumulant C2 associated with the graphs in Figs. [30 
are given by 



[< x 2 {t) >] = djt + d~ftT 1 - T u 



C 2 = T 1 - 2T 2 



(8) 



(9) 



The graphs in Figs. [30 and in Eqs. ((8][9|) appear in the 
same order from left to right. The identification of graphs 
with analytical expressions occurs according to existing 
rules. The effective vertex in these graphs is associated 
with the expression Go(V2(x 2 — £1), t 2 ~ ti)A(t 2 — As 
an example we give the analytical expression associated 
with Ti 

r i= J dt 2 ^ dh J d d x a J d d x b J d d x 2 J d d xi 

G (x a ,t- t 2 )G (x b ,t- i 2 )Go(v / 2 (x 2 - £1) ,i 2 - h) 
A(t 2 - OG^O, ii). (10) 



pt />t2 

T 1 = 2- d dt 2 <ftiA(t3-ti)(27T7ti)" 
Jo Jo 



■d/2 



(11) 



Similarly, the analytical expression associated with Ti x 
is obtained as 



Ti x = d~f2- d I dt 2 I dh[t-t 2 + 



10 Jo 
A(t 2 -ii)(2^i)"' i/2 

and consequently for the quantity 



2 4 



(12) 



H , -<hlT l = ,/ 7 2- d / dt 2 I dl } ( ■ '±1: 



10 Jo 
A{t 2 -t 1 ){2n 1 t 1 )- d ' 2 , 



(13) 



which yields the first-order correction to the mean-square 
displacement of the free end of DP in the skeleton expan- 
sion. 

The computation of Ti, d^tT\ 1 —T\ tX and C 2 in d = 1 
for large t results in 



y/2 Ao V 7T V t J 



+ ••., (14) 



7trx, -ria. 



try 



9 7 2 



2\/2A 2 

5 7 3/2 



16V2 4V2^A 



V~t- 



and 



V3^ A o V9 v^WA 



(15) 



"0 
(16) 

The computation of higher cumulants is similar, but is 
more time-consuming. The skeleton perturbation series 
permit to compute the quantities under consideration 
( [< x 2 {t) >] , C2, C3, • • • ) in the strong coupling regime 
at arbitrary dimensions. The skeleton perturbation ex- 
pansions possess, however, the deficiency which is due to 
the exponential dependence of the effective coupling con- 
stant on t. For sufficiently small time t the corrections 
are, however, small. While the connected graphs con- 
tain the factor exp(p c i), the disconnected ones the factor 
exp(kp c t) , where k is the number of connected parts of 
the disconnected graph. The expressions (|14H16p show 
that besides the exponential terms the expansion in the 
skeleton series occurs in inverse powers of A . Just above 
the transition (in d = 1 it corresponds to Ao ~ 0) the 
exponential terms are small, but the corrections to them 
are correspondingly large. 

We now will discuss the question how to proceed with 
the exponential terms, exp(p c t), in skeleton expansions. 
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A similar but however much simpler problem is the local- 
ization of a directed polymer in an external <5-potential 
U(x) = —u5(x) studied using the perturbation expansion 
in powers of u. The perturbation series of < x 2 (t) > is 
given in Fig. [5j 



of elimination of the exponential terms from the skeleton 
expansions, which physically is expected to correspond 
to the redefinition of the ground state, represents a new 
paradigm and remains unsolved yet. We explicitly com- 
pute the mean-square displacement of polymer end and 
the second cumulant of the free energy to the lowest or- 
ders of the skeleton expansions. 

A financial support from the Deutsche Forschungsge- 
meinschaft, the grant Ste 981/3-1, is gratefully acknowl- 
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FIG. 5: Perturbation series for DP in external potential. 

A straightforward summation of the series results in 
the exact expression of the mean-square displacement < 
x 2 (t) >= Z(t)/N(t), where the Laplace transforms of the 
numerator and denominator are respectively given by 

Zip) = ^(1-uG (0,p;0))-\ 

N(p) = -(l-uGofO.psO))" 1 . 
P 

The exponential terms in the perturbation expansion 
cancel in the exact expression. The DPRM can be con- 
sidered qualitatively as a localization on the randomly 
distributed ^-potentials, which suggests that an addi- 
tional partial summations of skeleton expansion might be 
necessary to eliminate the exponential terms, exp(kp c t) . 
However, the available terms in the skeleton expansions 
(I14j|16j) make it difficult to recognize the mechanism ac- 
cording to which such a resummation has to be carried 
out. Next-order terms in the skeleton expansions would 
simplify establishing the mechanism for performing such 
a resummation. While the bare perturbation expansions 
in powers of Ao are not alternating series, the skeleton se- 
ries (|14H16p are alternating ones. This property of skele- 
ton expansions supports indirectly the idea of a resum- 
mation. 

To conclude, we have performed partial summations 
of perturbation expansions of the directed polymer in 
disordered media to represent them as skeleton expan- 
sions in powers of the effective coupling constant, which 
corresponds to the binding state between two replicas 
or equivalently to the binding state of quantum particle 
in a delta potential. In the strong coupling phase the 
effective coupling constant exponentially increases with 
time. We argue that an additional resummation is nec- 
essary to cancel these exponential terms. The problem 
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